Theory of digital magneto resistance in ferromagnetic resonant tunneling diodes 
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We propose a ferromagnetic spintronic system, which consists of two serial connected resonant 
tunneling diodes. One diode is nonmagnetic whereas the other comprises a ferromagnetic emitter 
and quantum well. Using a selfconsistent coherent transport model we show that the current- 
voltage characteristic of the ferromagnetic diode can be strongly modulated by changing the relative 
orientation of the magnetizations in the emitter and quantum well, respectively. By a continuous 
change of the relative magnetization angle the total resistance exhibits a discrete jump realizing 
digital magneto resistance. The interplay between the emitter's Fermi energy level and the relative 
magnetization orientations allows to tailor the current voltage characteristics of the ferromagnetic 
diode from ohmic to negative differential resistance regime at low voltages. 
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I. INTRODUCTION 

The development of ferromagnetic dilute magnetic 
semiconductors (DMSs)^^^^ has opened the possi- 
bility of novel all semiconductor spintronic device con- 
cepts, in which the charge current can be modulated 
by the carriers sphii 7 - For example, spin dependent res- 
onant tunneling have been investigated both experimen- 
tally and theoretically in magnetic double barrier het- 
erostructures with either a ferromagnetic quantum well 
(QW), e.g., GaMnAs^^dLJ^i&iids or a paramagnetic 
QW, which exhibits a giant g-factor, e.g., ZnMnSe i 16 ' 17 ! 18 
By employing resonant interband tunneling an effective 
injection of spin-polarized electrons into nonmagnetic 
semiconductors has been demonstrated! 19 i 20 Moreover, 
enhanced tunneling magnetoresistance (TMR) has been 
predicted and found in double barrier magnetic tunneling 



junctions, 



9.14.15.21 



in which a nonmagnetic QW is sand- 



wiched between two magnetic electrodes. Recently, high 
magnetocurrents [relative current magnitudes for paral- 
lel (P) and antiparallel (AP) orientations of the magne- 
tizations] have been predicted in two coupled magnetic 
QWs^ 

Conventional nonmagnetic resonant tunneling diodes 
(RTDs) are technologically interesting due to their ex- 
treme high speed and low power performance. They al- 
low for novel circuit concepts based on their specific nega- 
tive differential resistance (NDR) behavior , 23 A logic gate 
named MOBILE (MOnostable-BIstable Transition Logic 
Element), which consists of two serial connected RTDs, 
a load and driver, has been proposed and realized by 
Maezawa and Mizutani i 24 i 25 The device is driven by an 
oscillating input voltage, which performs the transition 
between the mono- and bistable working point regimes. 
At low input voltages Kirchoff 's laws allow for only one 
stable dc working point. However, for high input voltages 
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two stable working points become possible due to the re- 
shaped current voltage (IV) charateristic of both the load 
and the driver RTD. Which of the two working points is 
actually realized depends on the difference of the load 
and driver peak currents. When the load peak current is 
higher than the driver one the working point voltage is 
high and vice versa. The whole device works actually as 
a comparator of the load and driver's peak current. For 
a detailed discussion of the operation principle see Refs. 

We have recently proposed^ that by replacing the 
driver by a magnetic RTD with a paramagnetic QW the 
circuit exhibits what we call digital magneto resistance 
(DMR): the output voltage jumps from low to high after 
the mono-to-bistable transition if the external magnetic 
field, which controls the Zeeman splitting in the QW, is 
higher than some threshold value. The threshold value 
of the magnetic field can be controlled by a gate voltage, 
which influences the peak current of the load. The pro- 
posed device performs a direct digital conversion of an 
analog magnetic signal and might be used as a fast mag- 
netic read head. DMR has actually been experimentally 
demonstrated in an earlier setup by shunting a metallic 
giant magnetoresistance (GMR) element to a nonmag- 
netic driver RTD, 27 having the advantage of being non- 
volatile upon the loss of power: the state of the device is 
stored in the magnetization direction of a particular fer- 
romagnetic layer of the GMR-element. Such nonvolatile 
devices are attractive for fast and reliable data storage, 
e.g, in random access memory applications 28 or for re- 
programmable logics, 29 in which the logical function of 
the circuit can be tuned by changing the magnetic state 
of the device. 

In this paper we lay down the physical principles of a 
nonvolatile ferromagnetic MOBILE. We propose to use 
a driver RTD, which comprises a ferromagnetic emitter 
and QW. By performing realistic selfconsistent calcula- 
tions of the IV-characteristics DMR is observed when the 
QW magnetization is tilted. For a proper choice of the 
emitter's Fermi energy the driver IV can be changed from 
ohmic to negative differential resistance behavior in the 
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low voltage regime depending on the relative magnetiza- 
tion orientation. 

The paper is organized as follows. A discussion of 
the system and the selfconsistent transport model is pre- 
sented in section [TXJ The simulation results are shown 
and discussed in detail in section HTTl and, finally, conclu- 
sions are given in section HVl 



II. MODEL 

The operation of the magnetic MOBILE is based on 
the change of the driver's peak current by applying an 
external magnetic field or by changing some magnetiza- 
tion direction in the device. In order to realize a non- 
volatile ferromagnetic MOBILE one can use the double 
barrier TMR structure^ in which a nonmagnetic QW 
is sandwiched between two ferromagnetic leads. For low 
voltages the structure yields large magnetocurrents, but 
for higher voltages, when the exchange splitting of the 
collector lead is shifted far below the band edges of the 
emitter lead by the applied voltage, the magnetocurrent 
becomes small. This means that for an exchange split- 
ting of the order of a few tens of meV the peak current 
is hardly influenced by the relative orientation of collec- 
tor's magnetization, since the peak voltage is already too 
high. Hence, hardly any DMR would be observed in such 
a TMR-device. 

Another setup to realize non- volatility would be to use 
a ferromagnetic QW instead of a paramagnetic one. In 
the proposed paramagnetic MOBILE^ the driver peak 
current decreases with increasing Zeeman splitting. In 
ferromagnetic materials it is easier to change the orienta- 
tion than the magnitude of the magnetization. If the ex- 
change splitting is strongly anisotropic a nonvolatile MO- 
BILE could indeed be realized by a ferromagnetic QW 
and by rotating its magnetization direction. However, 
in an isotropic case as considered here an additional fer- 
romagnetic lead is necessary to observe the DMR-cffcct. 
We consider a ferromagnetic emitter, since a ferromag- 
netic collector lead would again have little influence on 
the peak current. 

In order to investigate the characteristic physical ef- 
fects of such a structure as shown in Fig. [TJ we con- 
sider the emitter and QW to be made of a generic, 
moderately doped n-type ferromagnetic semiconductor. 
We also expect to observe similar effects in a p-type 
ferromagnetic semiconductor, since in the simplest ap- 
proach the heavy and light holes can be treated by an 
effective mass model^ analogous to our description of 
the conduction electrons here. In todays p-type ferro- 
magnetic DMSs, e.g., GaMnAs the Fermi energy can- 
not be chosen freely, since the ferromagnetic order ap- 
pears only at high hole densities^ However, for de- 
vice applications a decoupling of the ferromagnetic or- 
der and the doping would be advantageous. There have 
been several experimental reports on moderately doped 
n-type ferromagnetic semiconductors, e.g., HgCr2Se4, 32 




FIG. 1: (Color online) Left: The circuit configuration of the 
proposed ferromagnetic MOBILE. The load is a conventional 
RTD, whose peak current can be modified by an external 
gate voltage Vg- The driver device consists of a ferromag- 
netic RTD. The peak current of the driver is controlled by 
twisting the QW magnetization. Right: The schematic con- 
duction band profile of the ferromagnetic RTD [here made of 
a Ga(Al,Mn)N material system] used in the numerical simu- 
lations discussed in the text. 



CdCr 2 Se4r^ CdMnGeP 2 and most promising ZnO and 
GaMnNi^ i 35 ' 36 Many experimental results suggest room 
temperature (RT) ferromagnetism in transition metal 
doped GaN and ZnO. However, it is still controver- 
sial if the observed RT-ferromagnetism is intrinsic or 
due to some non-resolved precipitates. Several differ- 
ent mechanisms^ have been proposed theoretically to 
be responsible for the observed magnetic hysteresis. A 
thorough discussion of these issues is given in the review 
papers^£ and references herein. 

For the functioning of our proposed ferromagnetic MO- 
BILE we need a conduction band spin splitting of the or- 
der of tens of meV, regardless by which mechanism this 
exchange splitting is induced. In GaMnN it is generally 
believed that the exchange splitting of the conduction 
band is about 30-50 meVj 37 ' 38 In experiments the fer- 
romagnetic order sustains for dopings up to a few 10 18 
cm~ 3 , in thin layers of a few nm width.— It has also been 
shown theoretically that a small amount of anisotropic 
coupling in the 2D-Heisenberg model is sufficient to sta- 
bilize long range order at finite temperatures.— Two- 
barrier RTDs based on n-GaMnN have already been in- 
vestigated theoretically^ 8 - 

To be specific, we perform all numerical simu- 
lations for a GaMnN quantum well. The circuit 
diagram of the ferromagnetic MOBILE, the mate- 
rial composition, and the conduction band of the 
driver RTD is schematically shown in Fig. Q] Wc 
consider a two-barrier semiconductor heterostructurc, 
GaMnN/Ali^Ga^N/GaMnN/Ali^GazN/GaN, where 
a Ga concentration of about x — 17% is assumed in the 
barriers yielding a barrier height of about 200 meV. 40 
The QW is undoped, whereas the leads consist of 15 nm 
long n-doped layers, with n — 2.78 x 10 18 cm~ 3 in the 
magnetic emitter and n = 2.26 x 10 18 cm~ 3 in the collec- 
tor lead corresponding to a Fermi energy of Ef — 25 meV 
at the lattice temperature of T = 100 K. The magnetiza- 
tion M of the lead is considered to be fixed, whereas the 
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ferromagnetic QW is "soft" , which means that its magne- 
tization direction can be altered by an external magnetic 
field. 

In order to investigate the IV-characteristics of such 
a structure we follow the classic treatments of transport 
in nonmagnetic RTD ) 41 ' 42 ' 43 ' 44 ' 45 where coherent trans- 
port in the whole active device region (here the barriers 
and the QW) is assumed. Exploiting the symmetry of 
the Hamiltonian due to translations in the plane per- 
pendicular to the growth direction z of the heterostruc- 
ture, the Schrodinger equation can be reduced to a one- 
dimensional problem. By assigning the spin quantization 
axis to the fixed magnetization axis of the emitter lead 
the spinor scattering states ip z / E (z) with the spin quan- 
tum number a = ±1/2, (|, J.) regarding to the boundary 
condition of an incident plane wave form electron from 
lead i [= left (L), right (R)] with spin a' are then deter- 
mined in the effective mass envelope function approach 
by 



h 2 d l 



with 



2 dz m(z) dz 



U a {z) = E c {z) 



UAz) 



(z) = Etf (z), (1) 



e M (z). 



Here, m denotes the effective electron mass, and E is 
the total longitudinal energy of the electron (the sum 
of the potential and longitudinal kinetic energy), E c (z) 
denotes the intrinsic conduction band profile of the het- 
erostructure, e is the elementary charge, <f> the electro- 
static potential, A cx denotes the exchange splitting of 
the conduction band, a is the Pauli matrices vector, and 
e m is the unit vector of the magnetization. For a realis- 
tic simulations space charge effects have to be taken into 
account. The electrostatic potential 4> is obtained from 
the Poisson equation, 

±e(z)±<Kz) = -[n(z)-N d (z)], (2) 
dz dz eq 

where e denotes the static dielectric constant, eo is the 
permeability of the vacuum, Nj(z) is the fixed donor den- 
sity profile of the device, and n(z) is the electron density. 
The Poisson equation ([2]) has to be solved together with 
the Schrodinger equation (TTJ) in a selfconsistent way, since 
the quantum electron density is given by 



n(z) 



where 



-E 



hie) 



de Me) ^1 (z) 



1 



nh 2 



fc^Tln (1 + exp(/ii — e)) . 



(3) 



(4) 



Here, Uiy and v. 



denote the spin-dependent poten- 
tial energy and longitudinal group velocity of the elec- 
tron in the left and right lead, respectively, fii with 



/i/j = fiL ~ eVa is the chemical potential, where we as- 
sume that a voltage V a is applied to the right lead, and 
finally ks labels the Boltzmann constant and T the lat- 
tice temperature of the leads. Since we are interested in 
the generic properties of the proposed structure we ne- 
glect the effects of polarization charges at the interfaces, 
which appear in GaN due to spontaneous and piezoelec- 
tric polarization^ 

After obtaining the selfconsistent potential profile the 
current density is calculated in the framework of the 
Landauer-Buttiker formalism. By assuming parabolic 
bands and using the same effective mass for all layers of 
the heterostructure, the transmission matrix T a i a does 
not dependent on the transversal kinetic energy of the 
electrons. Hence, the current density j^J regarding to 
electrons which are incident with spin a from lead i and 
end up in lead j with spin a' can be obtained by gener- 
alizing the Tsu-Esaki formula 46 



+ j ernksT 
a = (2tt) 2 /i 3 



OO 

J deMe)TX>. (5) 



According to time-reversal symmetry T^J = T^. a , and 
the total current density is given by 



J 



^ ^ O'er' a 3 <t<t' ) ? 



(6) 



the difference of left and right flowing currents. 



III. NUMERICS AND SIMULATION RESULTS 

Following Ref. we numerically calculate the spinor 
scattering states in Eq. |T]) by applying the 4th-order 
Runge-Kutta (RK) scheme. Due to the forming of qua- 
sibound states in the QW the local density of states 
(LDOS) is strongly energy dependent. Therefore, we 
use an adaptive energy mesh for the numerical calcu- 
lation of the electron density by rewriting the quadra- 
ture, Eq. ([3]), into an initial value problem of an ordi- 
nary differential equation and solve it again by apply- 
ing a 4th-order RK scheme. In this way the computa- 
tional costs of the numerical integration are strongly re- 
duced (usually 400 energy points are necessary for a rel- 
ative accuracy of 10~ 4 , compared to 1200 grid points for 
an uniform energy mesh). To achieve fast convergence 
(usually 5-10 iterations steps) between the Schrodinger 
equations fT} and the Poisson equation we apply 
a predictor-corrector method^ The guess for the uni- 
formly discretized (Az = 0.1 nm) electrostatic poten- 
tial is obtained by a Newton-Raphson method, where the 
necessary Jacobian is estimated by using the semiclassical 
Thomas-Fermi approximation for the electron density4£ 
To save computational costs the particle density in the 
leads rii(z) is calculated semiclassically^i by 



rii(z) 



AV 
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1/2 



fj,i - U a {z)) 



(7) 
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FIG. 2: (Color online) Contour plots of the local density of 
states versus energy and growth direction z for parallel mag- 
netization alignment in the case of (a) equilibrium V a = 
and (b) at the peak voltage of V a = 0.12 V. The exchange 
splitting of spin up and down level is clearly visible for the 
quasibound ground state in the well. The solid lines indicate 
the self-consistent conduction band profile. 



where N c is the effective conduction band density of 
states and J- 1/2 is the Fermi-Dirac integral of order 
1/2. The transmission matrix T CTjCT / is obtained by the 
transfer-matrix technique. 41 Since the transmission func- 
tions are usually "spiky" we use an adaptive Gauss- 
Kronrod scheme for an efficient numerical calculation of 
the current density in Eq. (|5|). 

For all simulations we set the Fermi energy to Ef = 25 
meV (all energies are measured from the unsplit emitter 
conduction band edge), the lattice temperature to T = 
100 K, and we use the same exchange splitting in the 
QW and the emitter lead, A ex = 40 meV. We assume the 
spin up conduction band edge to lie energetically higher 
than the spin down one, which leads to a particle spin 
polarization of about -77.6 % in the emitter lead. An 
effective electron mass of m/mo = 0.228 40 (mo denotes 
the free electron mass) is used and the static dielectric 
constant is set to e = 9.54£ 

The contour plots of the local density of states (LDOS) 




AV i Voltage (V) 

FIG. 3: (Color online) Selfconsistent current-voltage charac- 
teristics of the magnetic driver- RTD for several relative orien- 
tations of the quantum well magnetization (indicated by the 
angle <p) at the temperature of T = 100 K. The solid black 
lines show the mirrored IV curve of the load- RTD for low and 
high input voltages, respectively; working points are indicated 
by circles. For these fixed low and high input voltages the out- 
put voltages are restricted to the intervals AV° ut ,i = 1, 2, 3. 
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FIG. 4: (Color online) Current-voltage characteristics of the 
magnetic driver-RTD in the non-selfconsistent case for differ- 
ent relative orientations tp of the quantum well magnetization 
at the lattice temperature T — 100 K. The peak voltage and 
currents are smaller when compared to the selfconsistent case 
of Fig. and there exists a common crossing points for all 
IV-curves. 



of the conduction electrons for P magnetizations align- 
ment in the considered heterostructure (see Fig. [T]) are 
shown in Fig. [DJa) in the equilibrium case (V a = 0) and 
(b) at the peak voltage of V a = 0.12 V. The forming of the 
exchange split spin up and down quasibound states is ap- 
parent. In the shown case of P alignment the spin up level 
lies energetically higher than the spin down one. The 
broadening of the first spin resolved quasibound states 
[at about 34 meV and 71 meV in Fig. Ufa)] in energy 
is much smaller than the exchange splitting, whereas for 
the next higher quasibound states (at about 164 meV 
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and 190 meV), which have a node in the middle, both 
spin levels already overlap. The solid lines indicate the 
selfconsistent conduction band profile. In the selfconsis- 
tent case parts of the applied voltage already drop before 
the first barrier but also beyond the second one as can 
be seen in Fig. Efb). This is in contrast to the non- 
selfconsistent case, where the applied voltage is assumed 
to drop linearly only in the active device region, i.e., in 
the QW and two barriers. 

The selfconsistent IV-characteristics for the struc- 
ture with the lateral dimensions indicated in Fig. [1] in 
the case of different relative orientations of the QW- 
magnetization (characterized by the angle ip) is shown 
in Fig. [3j For comparison, Fig. [4] displays the IV-curves 
for the non-selfconsistent case. Interestingly, in the lat- 
ter case there exists a single crossing point for all IV's 
with the same current 7* at some voltage V*. This fact 
suggests that the IV-characteristic for a specific angle ip 
might be written as a simple linear combination of the 
P and AP IV's: I V (V) = a(<p)I (V) + b(<p)I„(V) with 
a + b = 1. 

The crossing in the non-selfconsistent case follows from 
the linearity of the Schrodinger equation. For the pur- 
pose of the proof let us introduce a coordinate sys- 
tem, which we call the "0" -system, where the z-axis 
coincides with the magnetization direction of the QW- 
magnetization and the y-axis is given by the growth di- 
rection of the heterostructure. It is evident from the 
symmetry of the structure that instead of twisting the 
QW-magnetization by some angle — <p it is equivalent to 
assume a fixed magnetization direction in the QW and 
to rotate the emitter's magnetization by the angle <p. 
In the latter case only the boundary conditions for the 
Schrodinger equation (TTJ are changed when compared to 
the P alignment. The spinor wave functions in the emit- 
ter lead are of the plane-wave form, when they are repre- 
sented in a coordinate system (the 'V" -system) where the 
z-axis is given by the emitter's magnetization direction. 
The "ip" -system results from the "0" -system by rotating 
the latter around the y-axis by the angle ip. Hence, the 
spinor representations in the "0"- and "(^"-system are 
connected by {ip}o = D(—(p){ip} v , where 

n(,n\ - ( cos (v/ 2 ) sin(v3/2) \ . . 

°^>- (^-srm>/2) cos(^/2) )' (8) 

Due to the linearity of the Schrodinger equation (TTJ and 
since D(—ip) = cos((p/2)D(0) — sm(tp/2)D(ir) we can 
"divide" the problem into finding the solution for the 
two boundary conditions, {ip}o — cos(tp>/2){ r tp} v and 
Wlo = — sm(ip/2)D(Tr){^p} v , respectively. Represent- 
ing the latter in the 'V- coordinate system with its 
+ z-axis along the —z-axis of the "0" -system results in 
— sin(ip/2){ , ip} v . Thus, we can express the ampli- 
tude for the transmission of an incident electron with its 
spin aligned along the z-axis of the u (p" -system denoted 
by \Z(p, o) to a right moving plane wave state on the col- 
lector side with a spin eigenstate of the z-axis of the "0"- 
system \z ,a') as (zo,a'\z v ,a) = cos(tp / 2) {z , a'\z , a) + 
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FIG. 5: (Color online) Voltage-dependent magnetocurrent 
(MC) in the case of different relative orientations <p of the 
magnetization in the quantum well at the temperature of 
T = 100 K. 

sm(ip/2)(zQ,a'\z 7ri a). Since there is no spin precession 
for the P and AP case when the electron transmits 
the device region and since we assume that spin flip- 
ping scattering processes do not occur, the off-diagonal 
matrix elements of the amplitude matrix (zq, cr'\zo, a) 
vanish, whereas (zq, u'lz^, a) = for the diagonal ele- 
ments. With this the spin-dependent transmission func- 
tion, which is proportional to the squared transmission 
amplitude, can be written as 

T t',a a |(^0, cr'l^, cr)| 2 = 

/ cos 2 (f )\(z , T \z , T)| 2 sin 2 (f )\(z , 1 |)| 2 \ 
^skADKzoa l^,T>| 2 cos 2 (|)|(z ,l \z , l)| 2 J ' 

Using Eqs. |(SJ| and © finally yields the desired result, 

j v (V) - cos 2 (^/2)j (V) + sin 2 (^/2) j7r (V). (9) 

This relation is only valid in the non-selfconsistent case, 
since by including space charge effects, the selfconsis- 
tent electrostatic potential will become in general in- 
dependent and the above given considerations break 
down. We propose that the deviation from a common 
crossing point (/» , V* ) for all IV's might be used as a cri- 
terion of the relative importance of space charge effects, 
modulated by ip, in the device. As can be seen in Fig. [3] 
no common crossing point appears in the selfconsistent 
case. Since the applied voltage drops over a longer spa- 
tial region than in the non-selfconsistent case, as can be 
seen in Fig. [DJb), the selfconsistent peak voltages and 
currents are considerably higher. Hence, for our specific 
device setup space charge effects strongly influence the 
obtained IV-characteristics. 

The magnetocurrent (MC) for a particular magnetiza- 
tion orientation ip can be defined as follows: 

Mew = Mzjfr) , (10) 
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FIG. 6: (Color online) Current spin polarization as a function 
of the applied voltage for different orientations of the quan- 
tum well magnetization (indicated by angle <p) at the lattice 
temperature of T = 100 K. 



Our selfconsistcnt simulations reveal that the MC in- 
creases up to high values of about 800% (in the case of 
tp = 7r) for voltages right after the first NDR region of 
the AP IV-curve (V > 25 mV, see Fig. [J). Moreover, the 
current spin polarization with respect to the emitter spin 
quantization axis is given by 



■7'tt + in - hi - .in 

m +JU + JIT +311 



(ii) 



As shown in Fig. [5] Pj is strongly modulated when the 
QW magnetization is flipped from P to AP. At low volt- 
ages the polarization can be continuously changed from 
-93% for the P alignment up to +90% for the AP orien- 
tation. Hence, in the low voltage regime the ferromag- 
netic QW acts as a spin aligner. Interestingly, in the AP 
case the current polarization sharply slopes from +90% 
to -98% over the small first NDR voltage interval, which 
allows to use the device also as a voltage controlled spin 
switcher. 

To illustrate our results we introduce a model, which 
allows to give a qualitative estimate of the IV-curves for 
the extreme cases of parallel (P) if = and antiparal- 
lel (AP) ip = 7T magnetization orientation. In particular, 
the model allows for a better understanding for the in- 
fluence of the following three simulation parameters: (i) 
the energy difference for the lowest spin up and down 
quasi bound states 6 = — E± (it should be noted that 
S = 37 meV is a bit smaller than the exchange split- 
ting of the conduction band in the QW (A ex = 40 meV) 
due to the finite barrier height), (ii) the exchange split- 
ting A in the emitter lead, and (iii) the relative position 
of the Fermi energy level Ef in the emitter conduction 
band. By applying a voltage to the collector lead the 
quasibound states E a are shifted to lower energies. At 
first glance the voltage dependence of quasibound energy 
levels can be described by E a (V) — £^(0) — aV, where 
the parameter a is assumed to be voltage-independent 
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FIG. 7: (Color online) Illustration of the analytic model for 
the case of parallel (P) and antiparallel (AP) magnetization 
alignment. In the model the magnetization of the emitter is 
"soft" , whereas the QW-magnetization is assumed to be fixed. 
Part (a) shows the spin resolved conduction band profile in 
the emitter, the first quasibound states Ef and Ey in the 
quantum well, and the relative position of the Fermi energy 
level Ef. Part (b) displays the relative distances of the "on" 
and "off" -switching voltages of the spin- resolved components 
of the current on a 1/a voltage scale. The relative position 
between the "on" and "off" voltages (black filled circles) can 
be changed by moving the Fermi energy level. 



(a = 1/2 for the linear voltage drop and usually a < 1/2 
for the selfconsistent case) . For the following discussion it 
is more convenient to consider that the magnetization of 
the QW is fixed and, hence, in the AP case the magneti- 
zation of the emitter is flipped. This leads to completely 
identical IV-characteristics since in the AP case only the 
sign of the spin polarization of the current is changed but 
not its magnitude. The model is schematically illustrated 
in Fig. [3 

At low temperatures current can flow only if the qua- 
sibound states are dropped below the emitter Fermi en- 
ergy level, hence, the spin up and down currents are 
switched on at the voltages V° n — —(E^ — Ef). These 
spin-polarized currents are switched off, if the corre- 
sponding quasibound states are shifted below the emit- 
ter spin up and down conduction band edges, which 
leads to V^ oS = ±(E a - a A/2) for the P alignment 
and V££g = ±{E a + a A/2) for the AP case, respec- 
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tively. For higher temperatures the "switching on" re- 
lations will be thermally smeared on the order of a 
few ksT around the Fermi energy but the "switching 
off" relations are still valid, as long as inelastic scatter- 
ing can be neglected as assumed in our coherent trans- 
port model. From these relations it immediately follows 
that V^ s < y T P off < V[ oS < V T Ap ff for 6 < A and 
Vj^f < V p off < V p off < y T Ap ff for 6 > A, respectively. 
In the case of 6 < A the distances between these "off- 
switching" voltages are illustrated in Fig. [7]Jb) ; for 6 > A 
they are just given by interchanging (5 <-> A). 

First of all this model reveals that the peak voltage for 
the AP alignment is always higher than in the P case, 
since VJ^ S > V„ oS , as obtained in our simulated IV- 
curves Fig. [31 For the special case of 6 — A, which is 
approximately fulfilled in our simulations, the spin up 
and down current in the P case can be switched off at the 
same voltage. The relative position of the "on-switching" 
voltages on the voltage scale [indicated by the colored 
circles in Fig. Efb)] can be changed by the position of 
the Fermi energy level. By an appropriate choice of Ef, 
as indicated in Fig. [7^b), in the AP case the spin down 
current can be switched off before the spin up current is 
switched on. This leads to the NDR behavior in the low 
voltage regime, and the sharp current polarization drop. 
For the P alignment the total current is in the whole 
voltage range dominated by the spin down component. 
At low voltages up to 40 mV almost all current "flows" 
through the spin down channel but then also the spin up 
components starts to contribute to the total current, di- 
minishing the spin current polarization. Since both spin 
currents are "switched off" at almost the same voltage in 
the P case the IV-curve is nearly ohmic in the low voltage 
regime. This simple discussion shows that rotating the 
QW magnetization can drastically change the IV-curve. 

For the "intermediate" case of perpendicular orienta- 
tion (<p — 7r/2) (see Fig. [3]) two peak voltages appear 
with nearly the same values as those obtained for the 
P and AP alignment, respectively. However, most im- 
portant for observing DMR is that the peak current is 
already remarkably reduced when the relative magneti- 
zation orientation is tilted, say by an angle of a few tens 
of degrees, out of the P alignment. If we assume the 
load peak current is smaller than the driver peak current 
in the P case, the mono-to-bistable transition results in 
a low output voltage. However, if we tilt the magneti- 
zation orientation, the load and driver peak current be- 
come equal at some angle (p t h, which can be called the 
"threshold angle". For ip > y> t h a high output voltage 
is obtained in the bistable regime. Hence, the output 
voltage suddenly jumps from low to high after perform- 
ing the mono-to-bistable transition, effectively realizing 
DMR. The threshold angle can be controlled indirectly by 
properly tuning the load peak current, which is altered 
by an external gate voltage applied to the load device. 
Assuming a fixed low and high input voltage the output 
voltage is restricted to three different voltage intervals 
AV° U \ i = 1, 2, 3 as illustrated in Fig. H The high volt- 
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FIG. 8: (Color online) Blow up of the selfconsistent current- 
voltage characteristics displayed in Fig. [3] for low voltages. 
By replacing the load-RTD by a linear load resistance (indi- 
cated by the black solid line) allows to perform the mono-to- 
bistable transition by flipping the QW magnetization from P 
to AP. The working points are indicated by circles. The inset 
schematically shows the arrangement of the quasibound levels 
at the peak voltage in the AP case. 

age interval AVj™* is considerably separated from the low 
voltages intervals AV° ut and AV 2 out . This allows for a 
direct digital detection of the tilted QW magnetization. 

The question of how fast the mono-to-bistable transi- 
tion can be performed is closely connected to the sub- 
tle and still not clarified problem of determining the er- 
ror rates in MOBILES. The tilting of the QW mag- 
netization leads to a redistribution of the quasibound 
states, which usually takes place on the time scale of 
the order of a hundred of femtoseconds. The switch- 
ing time of RTDs, however, is limited by the "classi- 
cal" RC time constant, which is typically of the order 
of a few picoseconds^ In experiments the conventional 
MOBILEs randomly jump between high and low output 
voltage in the transition region and an erroneous transi- 
tion can occur due to parasitic capacitances or external 
electrical noised Transient studies of conventional MO- 
BILES based on an equivalent circuit model have shown 
that an error-free transitions with clock rise times on the 
order of the RC time of the RTD are possible if the out- 
put capacitance C ou t < (k — 1 ) Crtd r~ where Crtd is 
an average capacitance of the RTD and k is the ratio of 
load to driver peak current. Recently, conventional MO- 
BILES have been demonstrated to work up to frequen- 
cies as high as 100 GHz by employing a symmetric clock 
configuration for the input voltage. 52 This gives reason 
for a possible application of the proposed ferromagnetic 
MOBILE as a very fast "readout" of magnetically stored 
information. 

Fig. [8] shows a blow up of the selfconsistent IV- 
characteristics of Fig. [3] for low voltages. By rotating 
the QW magnetization the IV-curves change from ohmic 
to NDR-behavior. As illustrated in the inset of Fig. [5] 
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the spin down current becomes off resonance already be- 
fore the current can flow "through" the spin up quasi- 
bound state leading effectively to NDR. This interesting 
behavior can be used to perform the mono-to-bistable 
transition by twisting the QW magnetization instead of 
changing the input voltage from low to high as discussed 
above. For this, we assume to use a linear load resistance 
instead of the load-RTD in the circuit setup of Fig. [T] For 
an appropriately high input voltage two stable working 
points are obtained in the AP case, whereas only one 
crossing point appear in the load line diagram for the 
P alignment (see Fig. [5]) . With this again DMR can be 
realized as follows. Let us assume that at the beginning 
the circuit operates at the high voltage working point of 
the AP orientation. If we tilt the QW magnetization, 
suddenly at some threshold angle, the circuit is switched 
from the bistable to the monostable regime, leading to 
a discrete jump from high to low output voltage. This 
allows to detect "digitally" a disturbance of the AP mag- 
netization alignment. After this detection and after the 
recovering of the AP orientation the circuit will end up in 
the low voltage state of the bistable regime. By applying 
a small current pulse to the circuit it can again be re- 
set to the initial high voltage state. Another application 
could be a memory cell, in which the binary information 
is stored in the P and AP configuration and is read out 
by a small current pulse, which leads in the AP case to 
a voltage swing to the high voltage state. 

IV. CONCLUSIONS 

We have proposed a ferromagnetic MOBILE, where 
the driver-RTD comprises a ferromagnetic emitter and 
QW. By using a selfconsistent coherent transport model 
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